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We perform a simplified analysis of the edge excitations of the canted antiferromagnetic (CAF) 
phase of the v = quantum Hall state in both monolayer and bilayer graphene. Namely, we calcu- 
late, within the framework of quantum Hall ferromagnetism, the mean-field quasiparticle spectrum 
of the CAF phase neglecting the modification of the order parameter at the edge. We demon- 
strate that, at a fixed perpendicular component B± of the magnetic field, the gap A c d gc in the 
edge excitation spectrum gradually decreases upon increasing the parallel component By , as the 
CAF phase continuously transforms to the fully spin-polarized ferromagnetic (F) phase. The edge 
gap closes completely (A c d gc = 0) once the F phase, characterized by gapless counter-propagating 
edge excitations, is reached at some finite B^-dependent value By and remains closed upon further 
increase of By . This results in an gradual insulator-metal transition, in which the conductance 
G ~ (e 2 //i) exp(— Aodgo/r) grows exponentially with By in the range < By < By, while in the 
gapped CAF phase, and saturates to a metallic value G ~ e 2 /h in the F phase at By > By . This 
unique transport feature of the CAF phase provides a way to identify and distinguish it from other 
competing phases of the v — quantum Hall state in a tilted-field experiment. 

PACS numbers: 



I. INTRODUCTION 

Current transport experiments^^ provide compelling 
evidence for the interaction-induced nature of the ground 
states in monolayer (MLG) and bilayer (BLG) graphene. 
The most robust of them are observed in the quantum 
Hall regime at integer filling factors v corresponding to 
partially filled Landau levels (LLs). Most commonly, 
such states belong to the class of the so-called quantum 
Hall ferromagnets (QHFMs^pHSU _ bulk-incompressible 
states with spontaneously broken symmetry in the valley- 
spin space. 

The key physical challenge related to these states is 
identifying how exactly the symmetry is broken in a real 
system. One of the most intriguing questions concerns 
the nature the v — quantum Hall state with half-filled 
zero-energy LL, realized at the charge neutrality point. 
Its characteristic experimental signature in both MLG 
and BLG is the highly insulating behavior ^ 1 1 | 251 ^ of 
the two-terminal or Hall-bar longitudinal conductance, a 
strong indication that both bulk and edge charge excita- 
tions of the state are gapped. 

In Ref. a specific conclusion about the nature of the 
experimentally realized insulating v = state in BLG 
was made, namely, that it is a canted antiferromagnetic 
(CAF) phase of the v = QHFM, in which the spin po- 
larizations and Sx> of the valleys=sublattices=layers 
have equal projections on the direction of the total mag- 
netic field and are antiparallel in the perpendicular plane, 
see Fig. [I] and caption to it. This conclusion was based 
on the argument that CAF is the only phase on the 
generic phase diagram of the v = QHFM, obtained 
in Refs. 1221231 consistent with the transport data of 
Weitz et oZP and Velasco et a/P^ on dual-gated BLG de- 
vices, specifically, with the observation of the insulator- 



insulator phase transitions in the perpendicular electric 
field. 

Experimental verification of this conclusion requires 
concrete theoretical predictions for measurable quanti- 
ties, which would allow one to distinguish the CAF phase 
from other potential candidates. In this regard, the CAF 
phase of the v = state in both MLG and BLG is 
expected to exhibit a unique transport property in the 
tilted magnetic field. This property concerns the edge 
charge excitations of the CAF phase and was anticip ated 
in Ref. "231 based on the current understanding25ESH221 f 
the edge excitations of the antiferromagnetic (AF) and 
fully spin-polarized ferromagnetic (F) phases of the v = 
state and a general "by continuity" argument. 

Reiterating this argument, (i) the CAF phase contin- 
uously interpolates between the AF (9 S — n/2) and F 
(9 S = 0) phases, as the angle 29 s between the spin polar- 
izations sk,k' of the valleys=sublattices=layers is var- 
ied, skSk 1 = cos20 s , Fig. [TJ (ii) according to the ex- 
isting studies^^MlD the AF and F phases have gapped 
and gapless edge charge excitations, respectively (note 
that, at the same time, the bulk charge excitations of 
any phase of a generic QHFM are gapped); (iii) there- 
fore, by continuity, upon decreasing 9 S , as the CAF phase 
transforms to the F phase, the gap A ec i ge (9 s ) of the edge 
charge excitations of the CAF phase has to gradually de- 
crease and close completely once the F phase is reached, 
A odgo (0 S = 0) = 0. 

The optimal a ngle 29 s between the spin polariza- 
tions is controlle d 22 * 23 ^ by the ratio of the Zeeman en- 
ergy £z — I^bB, dependent on the total magnetic field 

B = Jb\ + Bp and the valley "isospin" anisotropy en- 
ergy uj_ = uj_(B±), dependent on the field component 
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FIG. 1: (Color online) Canted antiferromagnetic (CAF) phases of the v = quantum Hall state in monolayer (MLG) 
(left) and bilayer (BLG) (right) graphene. The spin polarizations sk and s K i of the valleys=sublattices in MLG and val- 
leys=sublattices=layers in BLG have equal projections on the direction of the total magnetic field (chosen as the z axis in the 
spin space) and are antiparallel in the perpendicular plane. In each valley, K and K' , the wave-functions of the zero-energy 
Landau level (e = LL) reside on either one of the sublattices, A or B, in MLG, and on either one of the sublattice, B or A, 
and therefore in either one of the layers in BLG. Thus for the v — state, the valley, sublattice, and layer degrees of freedom 
are equivalent: K o A, K' <H> B in MLG and K -o- B o (top layer), K' o A o (bottom layer) in BLG. Throughout the 
paper, we refer to the spin polarizations sk and Sj^ by their valley indices. 



B± perpendicular to the sample: 

ez 



cos 6>, 



2ujJ 



The anisotropy energy u±, defined in Sees. IV and 
VT[ originates from the Coulomb or electron-phonon in- 
teractions at the lattice scale, see Ref. [321 for details. 
As the key properties, u± is (i) linear in Bj_, uj_ ~ 
~ 1— 10B^[T]K, (a is some lattice spatial scale, 
yJc/(eB±) is the magnetic length, and we set 



e 2 a/l 



Ib 
h 



1 throughout the paper except for the c ondu ctance 
values e 2 /h), if the critical renormalizationj 22 * 32 * ^ are 
weak, or/and (ii) much greater than the Zeeman energy 
e z ± = HbB± for perpendicular field orientation, if renor- 
malizations are substantial. 

Therefore, in practice, the ratio ez/\u±\ can be effi- 
ciently changed and the transition from the CAF to F 
phase realized by tilting the magnetic field relative to the 
sample plane, Fig. [T] An unambiguous demonstration of 
the above predicted behavior of the edge gap A d gc of 
the CAF phase requires following the evolution of the 
system with varying the parallel field component B\\ at 
a fixed B±. This way the B^-dependent correlation en- 
ergies will not change in the process and the behavior of 
the CAF phase will be contrasted to that of the com- 
peting spin-singlet charge-density-wave or Kekule phases 
in MLG and fully layer-polarized or interlayer- coherent 
phases in BLG, whose bulk and edge gaps should not be 
sensitive (much, if at all) to the Zeeman effect. 

Thus, reexpressed in practical terms, at a fixed Bj_, 
the edge gap A e d ge of the CAF phase will gradually de- 
crease upon increasing By, as the CAF phase continu- 
ously transforms to the F phase. The edge gap will close 
completely at CAF-F phase transition point (which is 
of the second order at zero temperature) at some in- 
dependent finite value Bp determined from the condi- 
tion 



B 2 



Bjf 2 =2|ui(B]_)|. 



The edge gap will remain closed, A c d gc (By > By) = 0, 
upon further increase of By > B^ , as the system stays in 
the F phase. 

This behavior will manifest itself in the two-terminal 
or Hall-bar longitudinal conductance G as a gradual 
insulator-metal transition upon applying By . While in 
the CAF phase at < By < B* (e z < 2\u±\), the 
conductance should follow the Arrhenius activation law 
G oc (e 2 /h) exp(— Aedgc/7 1 ) determined by the edge gap 
Aedge (provided the contribution from the bulk with a 
larger gap Abulk > A d gc , nearly insensitive to By , is neg- 
ligible) and exhibit exponential sensitivity to By . Once 
the F phase is reached and upon further increase of By , 
at By > Bjj (ez > 2|uj_|), the conductance will saturate 

to metallic values G ~ e 2 /h due to conducting channels 
provid ed by the gapless counter-propagating edge exci- 
tationglSGSM] f the F phase. Ideally, if backscattering 
of the edge modes is negligible, the conductance in the 
F phase should be quantized as G — 2e 2 /h in MLG and 
G = 4e 2 /h in BLG, according to one and two channels 
per edge, respectively; in the two-terminal conductance, 
an extra factor of 2 arises from two edges, while in the 
Hall-bar longitudinal conductance it is due to the mode 
equilibration in the contacts. If partial backscattering 
is present, G will be lower and conductance fluctuations 
with varying By or other parameters can be expected. 

In this paper, we substantiate the above expectations 
by explicitly calculating the edge charge excitations of the 
CAF phase of the v = state in both ML G and BL G 
within a simplified approach. The majorit y! 15 * 28 * 30 ^ of 
the existing works on the edge excitations of the v = 
state in MLG neglect the modification of the bulk order 
parameter at the edge and calculate the mean-field quasi- 
particle spectrum. It was realized by Fertig and Brey^, 
on the other hand, who studied the edge excitations of 
the F phase in MLG, that the bulk order cannot be sus- 
tained at the edge due to the emergence of the finite ki- 
netic energy. They showed that, in fact, a "domain wall" 
is formed between the bulk and edge orders, where in 
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the latter electrons fully fill the hole branches of the edge 
spectrum. The proper lowest-energy edge charge exci- 
tations are then deformed configurations of the domain- 
wall texture in the valley-spin space that carry nonzero 
topological=electric charge. These excitatio ns are o f the 
same physical nature as the bulk skyrmion d 13 * 14 * 19 !, but 
have a crucially different energetics. 

As the comparison of the findings of Refs. 1151281301 
and Ref. [2^] shows, however, both approaches predict 
gapless edge excitations for the F phase. This suggests 
that, even though the former approach is not rigorous, 
its results would qualitatively agree with those of the 
latter for other phases, as well. For this reason and 
since the approach of Ref. [29] is technically consider- 
ably more sophisticated, in this paper, in order to pin- 
point the key physics, we follow the simplified approach 
of Refs. H5l28l30f3"Tl to study the edge excitations of the 
CAF phase of the v = state. Namely, we neglect the 
modification of the order parameter at the edge and cal- 
culate the mean-field quasiparticle excitations, and for a 
specific class of "armchair-like" boundaries. The analy- 
sis of the problem based on the generalization of the ap- 
proach^] of Fertig and Brey will be presented elsewhere^. 

We study the problem within the framework of 
QHFMisnp2H23] Before we proceed, we mention that an 
alternative to QHFMism approach to the v = state 
in MLG called "magnetic catalysis" was developed in 
Refs. |3"B144"3"1 While the approaches are arguably differ- 
ent, certain overlap between the results based on the two 
can be traced. In particular, the CAF phase of the v = 
state in MLG is predicted to undergo a similar evolu- 
tion in the tilted magnetic field within both the magnetic 
catalysis^! and QHFMisirP^l formalisms. 

The rest of the paper is organized as follows. In 
Sees. [TlfV] the mean-field excitations of the CAF phase 
of the v — state are studied for the case of MLG. In 

In 
In 



LL spacing, 



Sees Jill and III the model Hamiltonian is presented. 
Sec. |IV[ the mean-field bulk ground state is obtained. 
Sec. \Y\ the mean-field excitations of the CAF phase are 



obtained and their key properties are studied. In Sec. VI 



the findings of Sees. |H]|V] are generalized to the case of 
BLG. Concluding remarks are presented in Sec. |VII| 



II. HAMILTONIAN FOR THE ZERO-ENERGY 
LANDAU LEVEL IN MONOLAYER GRAPHENE 



Due to the formal equivalenc e 1 22 * 23 ! of the phase dia- 
grams for the v = QHFM in MLG and BLG, the results 
for the edge excitations of the CAF phase, as we show 
here, turns out physically the same as well. To keep the 
analysis clear, in Sees. [TlfV] we study the case of MLG 
in more detail and generalize the obtained results to the 
case of BLG in Sec. EH 

We start the analysis by writing down the projected 
Hamiltonian for the n = LL in MLG, valid at energies 
e <C v/Ib (v is the Dirac velocity) much smaller than the 



H — Hq + H[ Q + 



Hz 



#o = - E e (p)4 7 ^P' 



H; Q — 
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(2.1) 



(2.2) 



(2.3) 



h*=\ e 9a e ^r-^^M^y- 



a=x,y,z pi+ p> i=p2+ p> 2 



(2.4) 
(2.5) 



We consider a half-infinite sample occupying the x < 
half-plane and work in the basis of the bulk single-particle 
eigenstates 

\p\a) = \pX) <g> \a), A = K, K' and a =t, I, (2.6) 



\pK) = (W, 0,0,0), \ P K') = (0,0,0,V P ), 



3 i P y exp 



(x-x p ) 
211 



' Li, 



pi 



2 

B • 



(2.7) 



(2.8) 



of the n = LL in the Landau gauge A = (0,_Bj^a;,0). 
The components of the wave-functions \p\) are ordered 
as (^ka,iPkb,iPk'A,^K'b) in the KK' <E> AB valley- 
sublattice space. The states \p\cr) are characterized by a 
conserved one-dimensional momentum p along the edge 
and definite valley A and spin a quantum numbers. In 
each valley, K or K' , the wave- functions \pXa) reside on 
either one of the sublattices, A or B. 



For compactness, in Eqs. (2.2)-(2.5), we arrange the 
annihilation operators c v \ a of electrons in the states 
\p\o) into the spinors 

Cp = {CpKtiCpKlTCpK'tiCpK'lY 

in the direct product KK' s of the valley (KK') and 
spin (s) spaces; T a — t^ k <8> I s and S z — 1 KK <g) r| are 
the valley "isospin" and real spin operators, respectively, 
and : . . . : denote normal ordering of operators. 

Strictly speaking, the single-particle basis of the bulk 
states \p\cr) breaks down at the edge. As we show shortly 
in Sec. |III[ however, an important technical convenience 
is that, at energies e <C v/Ib of relevance to the low- 
energy theory of the v — QHFM, the effect of the edge 
can be incorporated perturbatively within the basis of 
the bulk eigenstates \p\a), and with minimal assump- 
tions about the edge properties. For a specific class of 
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FIG. 2: (Color online) The structure of the single-particle 
spectrum pertaining to the n = LL in a MLG sample with 
"armchair-like" boundary, neglecting the Zeeman effect. 



"armchair-like" boundaries we consider in this paper, this 
leads to the kinetic energy Hamiltonian Hq of the form 

As for the remaining terms in H, H lo describes the 
valley-symmetric screened Coulomb interactions, de- 



scribes the valley-asymmetric channels of the interac- 
tions, and Hz describes the Zeeman effect. The spin 
quantization axis z is chosen along the total magnetic 
field, which can have arbitrary orientation relative to the 
sample, see Fig. [T] 

The asymmetric interactions H lo arise from the ac- 
tual Coulomb interactions at the lattice scale or electron- 
phonon interactions with the optical phonon modes and 
may be taken as point in the real space. They are gener- 
ically characterized by two signed coupling constants 



g± — g x — g y and g Zl whose bare values scale as g, 
e 2 a and can undergo critical renormalizations^^ 



(0) 



32 36 



Explicitly breaking the valley symmetry, th e asymmetric 
interactions play a crucial role^MlH 22 * 23 * 44 * Ml j n selecting 
the favored ground state order of the v = QHFM. 
The standard expressions for the interaction matrix 

elements V p , lP , 2 and V PlP , 2 are provided in the Appendix. 

P1P2 P1P2 



III. "ARMCHAIR-LIKE" BOUNDARY 

In the effective low-energy Dirac theory for the elec- 
tron motion in MLG, valid at energies much smaller 
than the nearest-neighbor hopping amplitude t ~ 3eV 
between the carbon atoms, an edge of the sample is de- 
scribed by a boundary condition for the Dirac spinor 
V J = (4>KAjipKB,'4'K'A,ipK'B) t - As demonstrated in 
Ref. 1471 under general assumptions of preserved time- 
reversal and particle-hole symmetries, a generic edge of 
MLG is described by either an armchair or zigzag bound- 
ary conditions. 

In this paper, we restrict ourselves to the class of 
"armchair-like" boundaries, which, by definition, are de- 
scribed by an armchair boundary condition. Besides the 
actual armchair edge, more edge structures can belong to 
this class. A property of an armchair-like boundary, key 
to our consideration, is that it does not contain disper- 
sionless edge states in the absence of the orbital magnetic 



field. In the presence of the latter, this translates to the 
fact that the number of branches in the single-particle 
edge spectrum is equal to the discrete degeneracy of the 
bulk states \pXo-), i.e., four. 
Let 

\P~) = (lpKA,i>KB,i>K>A,'lpK'B) 

be an ex act solutio n of the boundary problem in a mag- 
netic fiel d 28 * 30 * 48 * 49 ! for the Dirac equation per given spin 
projection a, characterized by a negative energy — e(p) < 
(hole branch) and pertaining to the n = LL. The 
kinetic energy e(p) — > is flat in the bulk (p < 0) and 
grows a t the edge (p > 0); its exact functional, although 
obtaine d 28 * 30 * 48 * 4 ^ , is not essential right now. The eigen- 
state 



KA, 



'IpKB^K'A, -IpK'B)- 



(3.1) 



with the positive energy e(p) > (particle branch) is 
obtained from \p—) by the particle-hole transformation, 



-e : ip. 



XB 



XB; 



A = K,K' 



(3.2) 



By the assumption of an armchair-like boundary, there 
are no other single-particle states besides |p±) at energies 
|e| -C v/Ib- The second-quantized Hamiltonian for the 
kinetic energy reads 



H o=Yl e (p)( c l+<y c p+° - c l-a c P-<y)> 



(3.3) 



where c p ± a are the electron annihilation operators for 
the states |f>±c) = |p±) S) |<r), with spin included. The 
spectrum is shown in Fig. [2j 

Deep in the bulk, where the kinetic energy e(p) — > 
vanishes, the exact eigenstates |p±) must evolve into the 
linear combinations of the bulk states \pK) and \pK') 
[Eqs. (2.7) and (2.8)] that are related by the particle-hole 



symmetry, 
|P±) 



V2 



(\pK) T \pK')) at e(p) -> 0, 



(3.4) 



The electron operators are related accordingly, 

c P ±<r = —7^(cpKa T CpK'o-) at e(p) -> 0. (3.5) 
v2 

Immediately at the edge, where e(p) ~ v/Ib, the bulk 
states \p\cr) [Eqs. (2.6)-(2.8)] are not well-defined. How- 



ever, not too close to the edge, where e(p) <C v/Ib, one 
may treat the kinetic energy as a perturbation, neglect- 
ing the modification of the wave- functions but taking into 
account the energy splitting of the states \p\cr). In this 
case, one may still use the relation ( |3.5[ ). Substituting 
Eq. (3.5) into Eq. (3.3), we obtain an approximate ex- 



pression (2.2) for the kinetic energy in the basis of the 
bulk states |pAcr), valid not too close to the edge, so that 
t{p) < v/l B . 
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Thus, at energies below the LL spacing v/Ib, the ef- 
fect of an armchair-like edge on the n = LL states 
can be taken into account perturbatively while remaining 
within the basis of the bulk eigenstates |pAcr). The ef- 
fect amounts to an effective "Zeeman" field — e(p)(l, 0, 0) 
along the x direction in the K /T'-isospin space that hy- 
bridizes the \pKa) and \pK'a) states, favoring the occu- 
pation of the states \p— a) of the hole branch, Eqs. (3.1 1, 
(j3.4h , and (3.5). Note that the direction of this "Zeeman" 



field in the xy isospin plane could be chosen arbitrary 
due to the freedom of choice of the phase factor in the 
particle-hole transformation (3.2). 



IV. BULK GROUND STATE 

The strongly interact ing v = state is described by 
the theory of QHFMisnP^HSIl. I n this Section, we briefly 
recover the results of Refs. 1221231 for the bulk ground 
state, pertaining to the CAF phase. 

One constructs a Slater-determinant state 



P V Act / 



E< AV IX0> C pA'ct< ) |0>, 
\\'o' 



(4.1) 

in which two electrons per each orbital p of the n = LL 
occupy arbitrary mutually orthogonal states Xa,b in the 
KK' (g) s space. It is straightforward to show that <J/ is an 
exact eigenstate of the Hamiltonian (2.3) of the SU(4)- 
symmetric interactions for any choice of the spinors Xa.bj 

For a wide class of repulsive interactions, one can ex- 
pect the eigenstates ^ to be exact ground states by the 
Hund's rule argument. This is the main assumption of 
the QHFMism theory, also employed in this paper. 

Thus at the level of symmetric interactions, the ground 
state is known exactly, but it is highly degenerate. This 
degeneracy is uniquely parameterized by the order pa- 
rameter matrix 

P = Xaxl+Xbxi (4.2) 
which satisfies the properties of a projection operator, 
pt = P p2 = p tT p = 2 ^ ( 4 3 ) 

The favored order P is determined by the effects that 
explicitly breaks the SU(4) symmetr y in the KK' ® s 
space: valley-asymmetric interactions (2.4) and the Zee- 



man effect ( 2.5 ). These effects are taken into account per- 



turbatively by calculating their energy expectation val- 
ues, 

£{P) = £ (P)+£ Z (P), (4.4) 
£ (P) = (*|flio|*)/JV = 

= J 52 u a {tv 2 [T a P] -tr[T a PT Q P]},(4.5) 



Here N = J2 P 1 i s the number of orbital states, equal to 
the number of flux quanta threading the sample. The 
asymmetric interactions result in the isospin anisotropy 
energy £ (P), characterized by two signed energies 



Uj_ 



ff-L 

2nl 2 B 



9z 



2tt/|' 



(4.7) 



Minimization of the energy £(P) of the SU(4)- 
symmetry-breaking effects for arbitrary values of uj_, z 
and tZi resulting in the generic phase diagram for the 
v = QHFM in MLG and BLG, was carried out in 
Refs. 1221231 In this paper, we will be interested in the 
canted antiferromagnetic (CAF) phase, argued in Ref.l23l 
to be realized in the insulating v = state of the real 
BLG. 

The CAF phase is realized when the isospin anisotropy 
£o(P) alone favors the AF phase. This occurs for u z > 
—uj_ > 0, the condition assumed to be satisfied in the 
rest of paper. In this case, in the presence of the Zeeman 
effect, the energ y £(P) is minimized by either the CAF 
or F phase^HM with \a — \K) <g> |sj<-) Xb 
in which the spin polarizations 



\K') 



\SK' 



sk,k' = (i sin 9 S cos (p s> zL sin 8 S sin tp s , cos 6 S ) (4.8) 

of the valleys=sublattices in MLG and val- 
leys=sublattices=layers in BLG have equals projections 
on the direction of the total magnetic field and are 
antiparallel in the perpendicular plane, Fig. [T] 
The energy 

£ (P) = —u z — u± — u± cos 28 s — 2ez cos 6 S 

of this family of states is minimized at 



\ 1, e z >2|u x |(F). 



(4.9) 



I.e., for ez < 2\u±\ the ground state is a CAF phase with 
the optimal angle 28 s between the spins and energy 



£ 



CAF 



2«x 



while for ez > 2\u±\ it is a fully spin-polarized F phase 
with the energy 

£ F = 2\u ± \-u 2 . 

The CAF phase has a U(l)-degeneracy according to 
the choice of the spin orientation ip s in the plane per- 
pendicular to the total magnetic field. For the calcula- 
tions below, we will assume a specific orientation ip s = 0, 
in which case the order parameter (4.2) of the CAF/F 
phases takes the form 

P= ^1 ® (1 + cos 9 s t z ) + ^t z <E) sin 9 s t x (4.10) 
As discussed in the Introduction, the obtained bulk 



£z{P) = {*\H z \$>)/N = -e z tr[S z P] 



(4.6) order parameter (4.10) cannot be sustained at the edge 
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due to the emergence of the finite kinetic energy and will 
be necessarily modified at p, such that e(p) > £(P). For 
the reason mentioned there, however, in this paper, we 
will neglect this fact and use the bulk order parameter P 
at all momenta p, even where the kinetic energy is not 
negligible anymore. 



V. MEAN-FIELD EXCITATIONS OF THE 
CANTED ANTIFERROMAGNETIC PHASE IN 
MONOLAYER GRAPHENE 

Mean-field excitations are obtained by performing de- 
coupling of the interactions in Eqs. (|2.3[) and (|2.4|), 



#i,mf = ~ y^&l Ac p , 



(5.1) 



A = u a P- J2 u a (T a tr{PTa}-T a PT a ), (5.2) 



where u — J2 p , V pp f is the exchange energy of the sym- 



metric interactions (2.3) (we discard the P-independent 
Hartree energy of the symmetric interactions) and the 
isospin anisotropy energies uj_ z of the asymmetric inter- 
actions (2.4) were defined in Eq. (4.7). 



Inserting Eq. (4.10) into Eq. (5.2) and discarding the 



trivial terms tx 1(8 1, we obtain for the mean-field poten- 
tial of the AF/CAF/F phases 

A = A *i ® t z + A zxTz ® t x , (5.3) 

A 0z = -(u Q +u z -2\u ± \)cos6 s , (5.4) 

A zx = -(u Q +u z + 2\u_ L \)sind s . (5.5) 

The full mean-field Hamiltonian takes the form 

Hml = H o + -ffi, m f + Hz = ^ cphcp, 



h = -e(p)r x <E)i-(e z + A 0z )i<E)T z 



>t x . (5.6) 



The single-particle Hamiltonian h is straightforwardly 
diagonalized. This yields the four branches 

E±±{ P ) = ±VHp)±(ez + A 0z W + Al x (5.7) 

of the mean-field spectrum of the AF / CAF /F phases of 
the u = state in MLG, plotted in Fig. [3] 

At the edge, the branches (p) and E^ + (p) = 

— (jp) come closest to each other and the gap A e d ge in 

the edge excitation spectrum is determined by the branch 
minimum, A cdgc = 2£ , + _(p ), where d£ + _ (p)/dp\ Po = 
0. We obtain 



^V 1 -^) 2 ' ^<2KI(CAF), 
0, e z >2\u x \(F), 



where 



A AF = u + u z + 2\u ± \ 



(5.9) 



is the edge and bulk gap of the AF phase at ez = 0. 

The bulk gap Abulk = 2E + ±(p — > —oo) is obtained by 
setting e(p) = in Eq. (|5.7[) and equals 



A bulk = 2^/(A 0z + e z y+Al x = 



A AF , 
a af 



-2{e z 



ez < 2\uj_\ (CAF), 
■2\u±\), e z >2\u x \(F). [oAU) 



Equations &5.7L (5.8), and (5.10) for the spectrum of 



the mean- field excitations of the AF/CAF/F phases in 
MLG constitute the main result of the work. Below we 
discuss their key properties and demonstrate that they 
do lead to behavior anticipated in the Introduction. 



The mean-field potential (5.3) of the CAF phase is 
a mixture of the ferromagnetic (A 0z ) and antiferromag- 
netic (A zx ) components; their relative value is fully con- 
trolled by the ratio ez/\u±\, which determines the angle 
29 s between the spins polarizations sk,k' [Eq. (4.9)]. 



A 



Cclff 



2A ? 



As discussed in the Introduction, efficiently changing 
the ratio ez/\u±\ in the experiment requires tilting the 
magnetic field, Fig. [l] Let us discuss the evolution, plot- 
ted in Figs. [3] and [4j of the edge A ec ige [Eq. ( 5.8 )] and bulk 
Abulk [Eq. p,10| )]""gaps upon changing the Zeeman en- 
ergy ez at fixed values of the interaction energies uo t ± tZ . 
Practically, this corresponds to fixing B± and applying 

According to Eq. (5.8), the edge gap A d go is equal 
to twice the AF component A zx of the mean-field po- 
tential. In the theoretical limit of vanishing Zeeman en- 
ergy ez — 0, the phase is purely AF; the F component 
Aqz = is absent and the edge and bulk gaps are equal, 
A cdgc = Abulk = A AF , Fig. |3[a). Upon increasing e z , 
the edge gap A d g c = 2A ZX of the CAF phase gradu- 
ally decreases as sin S = yjl — (^fj) 2 : as tne an gle 20 s 

between the spins decreases and the CAF phase continu- 
ously transforms to the F phase, Fig. [3jb) . The F phase 
is reached at ez — 2\u±\ and persists upon further in- 
crease of ez > 2|u_i_|; the AF component turns zero and 
the edge gap closes, A c d g c = A zx = 0, Fig.^c). The F 
phase is characterized by the gapless counter-propagating 
edge excitations with opposite spin projections and spec- 
tra ±[e(p) — ez], in accord with earlier findings 28 " 30 . 

At the same time, according to Eq. (5.10), the mean- 
field bulk gap Abulk = A AF of the CAF phase, e z < 
2|uj_|, does not d epen d on the Zeeman energy ez and is 
equal to the gap (5.9) of the AF phase. In the F phase, 
ez > 2|ujJ, Abulk grows with ez due to the Zeeman 
effect. 

These findings explicitly confirm the expectation of 
Ref. [531 for the properties of the edge excitations of the 
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FIG. 3: (Color online) The mean-field spectrum of the (a) antiferromagnetic (AF), (b) canted antiferromagnetic (CAF), and 
(c) ferromagnetic (F) phases of the v = state of a monolayer graphene (MLG) sample with armchair-like boundary (see 
Sec. |III| for a definition). The spectrum consists of four branches E±±(p) [Eq. ( |5.7| l]. The edge gap A c d ge is determined by the 

shortest distance between the E^ (p) and E- + (p) — — E+-(p) branches. The edge gap A c d gc is maximal in the AF phase (a), 

has a smaller value in the CAF phase (b), and vanishes in the F phase (c), wh ere g apless counter-propaga ting e dge excitations 
with opposite spin projections emerge. The evolution of the edge A c d gc [Eq. (5.8l] and bulk Abuik [Eq. ( 5. 10 1] gaps with the 
Zeeman energy ez at fixed interaction energies «o,±,z is presented in Fig. [4] 




2|«J 



AF 

i r 



CAF 



T I / / X 



2|«J 



FIG. 4: (Color online) (Top) The evolution of the edge A e d ge 



[Eq. ( 5.8 1] and bulk Abuik [Eq. (5.10l] gaps in the mean-field 
excitation spectrum [Eqs. ( |5.7[ ) and ( 6.14[ ) and Figs. [3] and [6] 
of the AF/CAF/F phases of the v = state in both MLG 
and BLG upon changing the Zeeman energy ez at fixed val- 
ues of the interaction energies Mo,J_,z- In practice, such evo- 
lution is realized by varying a parallel field component B\\ of 
the magnetic field at a fixed perpendicular component B±, 
Fig. [I] The edge gap A c d gc , determin ed b y the AF com- 
ponent A zx of the mean-field potential J5T3J) , is maximal in 
the AF phase at ez = [Figs. [3|a) and [6ja)], gradually de- 
creases upon increasing ez in the range < ez < 2|w_l|, while 
in the CAF phase [Figs. |3|b) and[6|b)], and vanishes at the 
CAF-F phase transition at ez — 2\u±\ and upon further in- 
crease of ez > 2 1 iix I, as the system stays in the F phase 
[Fig. [3jc) and Fig. [6jc)]. The bulk gap Abuik is constant in 
the AF/CAF phases at ez < 2\u±\ and grow in the F phase 
due to the Zeeman effect. (Bottom) The corresponding evo- 
lution of the spin polarizations s K K > [Eqs. (4.8 1 and (4.9 1] of 
the valleys=sublattices=layers in the AF/CAF/F phases. 



CAF phase of the v — quantum Hall stated The re- 
sulting physical behavior and experimental implications 
were discussed in the Introduction. 



VI. MEAN-FIELD EXCITATIONS OF THE 
CANTED ANTIFERROMAGNETIC PHASE IN 
BILAYER GRAPHENE 



The findings of Sees. P|V| for MLG are straightfor- 
wardly generalized to the case of BLG. The key extra 
feature in BLG is that both n = and n — 1 mag- 
netic oscillator states belong the e = LlP^, resulting 
in its additional two-fold degeneracy. The bulk eigen- 
states \pn\a) of the e = LL in the Landau gauge are 
characterized by momentum p, valley A = K,K', spin 
a =~f, 4,, and n = 0, 1 quantum numbers. Each orbital p 
is therefore eightfold-degenerate. 

The projected Hamiltonian for the e = LL in BLG, 
valid at energies e <§; l/{ml%) below the LL spacing (m 
is the effective mass of the quadratic spectrum), has the 
form 



H = H a 



H< 



Hz, 



Ho = - y^^njp^ynTxCpn, 



(6.1) 



(6.2) 



Pl+P' 1 =P 2 +P 2 

a=Q,x,y,z Pi+p' 1 =p 2 +p' 2 
X : [^p 1 n 1 TaCp 2 n 2 ][Cp' in ' i TaCp 2 n 2 ] '■ 



(6.3) 



n 2 p 2 
i-' 2 P' 2 



Hz ^ Cpri^ zCpnj 



(6.4) 



pn 



IpO+O-)/ |/>l+<r> 




FIG. 5: (Color online) The structure of the single-particle 
spectrum pertaining to the e = LL in a BLG sample with 
"armchair- like" boundary, neglecting the Zeeman effect. 



Similarly to the considerations of Sec. |III| at energies 
e <C l/(mZ|j), the effect of the edge of a BLG sample 
can be incorporated perturbatively within the basis of 
the bulk eigenstates \pn\o) of the e = LL. For BLG, 
we also restrict ourselves to the case of "armchair-like" 
boundaries, which, by definition, are described by an 
armchair boundary condition in the two-band model of 
BLG^ and, consequently, do not contain dispersionless 
edge states in the absence of the magnetic field. In this 
case, the LL spectrum has 8 branches pertaining to the 
e = LL. Deep in the bulk, where the kinetic energy 
vanishes, the eight eigenstates \pn±a) of the boundary 
problem must evolve into the linear combinations of the 
bulk states \pn\a) that are related by the particle-hole 
symmetry, 



\pn±a) 



1 



(\pnKcr) T \pnK' a)) at e n (p) -> 0, (6.5) 



V2 



(CpnKcr T c pnK'a) a t e 7i(p) 



0. 



(6.6) 



This consideration simultaneously proves that only the 
different valley states \pnKa) and \pnK'a) with the same 
n and a quantum numbers get hybridized by the edge. 
Thus, the groups \p0\a) and [pi Act) of the bulk states 
evolve into two uncoupled sets of branches with kinetic 
energies ±eo(p) and ±ei (p), as shown in Fig. [5] which 
also justifies labeling the exact eigenstates \pn±a) by the 
n = 0, 1 quantum numbers. Since the n = 1 states have 
a larger spatial extent than the n = states, the energies 
(arising from hybridization) are related as ei(p) < eo(p)- 
The exact kinetic energy Hamiltonian for the states per- 
taining to the e = LL reads 



pn-\-(T 



pn 



Substituting Eq. (6.5 1 into Eq. (6.7), we arrive at the 



approximate expression (6.2) for kinetic energy Hamilto- 



nian within the basis of the bulk eigenstates \pnXa) . 

The terms (6.3) and (6.4) describe the valley- 
symmetric and asymmetric interactions, respectively. 
Due to the difference in the orbital wave-functions of the 
n = and n = 1 states, both interactions are necessarily 
anisotropic in the 01-subspace. The valley-asymmetric 
interactions (6.4) are point to a good approximation. 



To keep the analysis simpler, we consider the valley- 



symmetric interactions (6.3) as point, as well, in which 
case uq = <7o/(ti^b) m the equations of Sec.[Vl Consider- 
ing finite-range symmetric interactions would somewhat 
modify the spectrum quantitatively, but not qualita- 
tively. The expression for the matrix element V n , lP , 1 ' n ? p , 2 
is provided in the Appendix. 

The technical steps leading to the edge excitations are 
analogous to the case of MLG. One first obtains the bulk 
ground state of the v = state within the framework of 
QHFMism. As demonstrated in Refs. USlT^l the asymme- 
try of the interactions in the 01-subspace favors electrons 
to arrange into 01-singlet pairs with identical content in 
the KK' <E> s space within a pair. Precisely, the states 



* = II Ew<U £<AV| X6 >ct 



pn\' <r' 



|o>, 



deliver the variational minimum among the QHFM states 
in the full 01 <g> KK' (g> s space for the valley-symmetric 
interactions H 10 . It is also possible to show that (6.8) is, 
in fact, an eigenstate, 

H io V = S *. 

However, the mutually orthogonal spinors Xa.b in the 
KK' (X> s space can still be arbitrary . Th us, at the level 
of the valley-symmetric interactions (6.3), the (assumed) 
ground state * is degenerate. 

To find the preferable order P = XaXL + XbXb m 
the KK' <S> s space, one calculates the energy of the 
symmetry-breaking terms, the asymmetric interactions 
(6.4) and Zeeman effect (6.4), 



£(P) = £«{P) + £ Z {P), 
£o(P) = (tt|fli o |*>/(2A0, 
£z{P) = {^\H Z \^)/{2N). 



(6.9) 
(6.10) 
(6.11) 



With the definitions 



Uj_ = U X = Uy — 



9± 



9z 
7:11 



- a ). (6.7) the functional dependencies of the isospin anisotropy 
£o{P) and Zeeman £z(P\ energies on P app ear t o be 
identical to those in MLGpS, as given by Eqs. (4.5) and 



As in MLG, not too close to the edge, wher e e n (p) <C (4.6). Consequently, the phase diagram for the v = 



\/{ml B ), one may treat the kinetic energy (6.7) as a 



perturbation and still use the relations (6.5) and (6.6) 



state in BLG, obtained by minim izing the energy £(P), 
turns out formally identical 22 * 23 ^ to that in MLG. 
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(a) AF (e z =0) ' 






t 

Abulk=A e dge 




* 


p 








FIG. 6: (Color online) The mean-field spectrum of the (a) antiferromagnetic (AF), (b) canted antiferromagnetic (CAF), and 
(c) ferromagnetic (F) phases of the v — state of a bilayer graphene (BLG) sample with armchair-like boundary. The spectrum 
(6.14 1 consists of eight branches E n ±± (p) [Eq. ( 6.14 1] - two uncoupled sets of four branches pertaining to the n — (blue) and 

and bulk Abuik [Eq. ( |5.10| l] gaps are the same for n — and n = 1 sets and 



(green) states. The edge A d gc [Eq. (|5.8 



the same as in MLG; their evolution with the Zeeman energy ez at a fixed anisotropy energy u± is presented in Fig. [4] 



Owing to this correspondence, the results for the exci- 
tations of the AF/CAF/F phases are basically the same, 
as well. The mean- field Hamiltonian equals 



Hmf = H + -ffi.mf + Hz 



pn 



-e n (p)T x <S>i - {e z + A 0z )l( 



-A, 



(6.12) 



(6.13) 



Diagonalizing h n for each n = 0,1, we obtain the eight 
branches 

E n±± (p) = ±y/[e n (p)±(e z + A 0z )] 2 + A* x , (6.14) 

of the mean- field spectrum of the AF / CAF /F phases of 
the v = state in BLG, shown in Fig. [6] The condi- 
tion |j49j) for the optimal angle 6 S of the AF/CAF/F 
phases, the expressions for the bulk order parameter 
P [Eq. (|410j)1, components A 0z [Eq. Kty] and A zx 



Eq. 5.5)] of the mean-field potential A are the same as 



in MLG. 

The only essential difference of BLG from MLG is that 
due to the extra two-fold degeneracy of the e = LL in 
BLG, the number of branches is doubled: there are two 
uncoupled sets, n — and n = 1, of four branches and 



for each set the spectrum (6.14) is the same as in MLG, 
compare with Eq. ( 5.7 1 . The n — and n = 1 sets are not 



coupled by the inter actio ns because the order parameter 
l 01 (g) P of the state (6.8) is singlet in the 01-subs pace. 
The edge A odgo [Eq. (|5.8|)] and bulk A bu ik [Eq. (pUO)] 



gaps of the AF/CAF/F phases in BLG are the same as in 
MLG and the same for n = and n = 1 sets, Fig. [6] Also 
identical is the evolution of the A cc j gc and bulk Abuik g a ps 
with varying ez /\u±\ by tilting the magnetic field, Fig. [4] 
Since in the F phase (ez > 2|u_l|), there are two gapless 
edge channels for each direction, n = and n — 1, for 
negligible backscattering, the two-terminal or Hall-bar 
longitudinal conductance G BLG = 4e 2 /h of the F phase 
is quantized at twice the MLG values. The resulting 
physical behavior and experimental implications are thus 
essentially identical to those in MLG and were discussed 
in the Introduction. 



VII. CONCLUSIONS 



In conclusion, we studied the charge excitations of 
the canted antiferromagnetic phase of the v — quan- 
tum Hall state in monolayer and bilayer graphene within 
a simplified approach, namely, neglecting the modifica- 
tion of the order parameter at the edge and calculat- 
ing the mean-field spectrum for a class of "armchair- 
like" boundaries. We explicitly demonstrated that the 
gap in the edge spectrum of the canted antiferromag- 
netic phase monotonically decreases upon decreasing 
the angle between the spin polarizations of the val- 
leys=sublattices=layers, as the canted antiferromagnetic 
phase continuously transforms to the fully spin-polarized 
phase, Fig. [l] |3j |4j and [6j In practice such evolution 
is realized by tilting the magnetic field and, in order 
to contrast this behavior to that of the competing spin- 
singlet phases, the evolution with varying Bii and fixed 
B± should be traced. The edge gap closes completely 
as the fully spin-polarized phas e is reached and gapless 
counter-propagating modes^^SHini emerge. This results 
in an gradual insulator-metal transition, where the con- 
ductance grows exponentially in the insulating canted an- 
tiferromagnetic phase and saturates to metallic values in 
the fully spin-polarized phase. 



These unique edge transport properties of the canted 
antiferromagnetic phase, earlier anticipated in Ref. 1231 
provide a straightforward way to identify and distinguish 
it from other competing phases of the v — state in the 
experiment. The data of the recent tilted- field experi- 
ment on bilayer graphene by Maher el a/P^l - exponen- 
tial growth of the longitudinal Hall-bar conductance G xx 
over a finite range of applied Bn at fixed B±, followed by 
the saturation to metallic values G xx ~ e 2 /h - are fully 
consistent with this predicted behavior, lending crucial 
support to the conclusion of Ref. [23] that the insulat- 
ing v = q uantum Hall state realized in real bilayer 
graphen e 1 8 * 1 ^ is canted antiferromagnetic. 
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Appendix: Interaction matrix elements 



The interaction matrix elements in Eqs. (2.3), (2.4 1, 



(6.3), and (6.4) are given by the standard expressions 



yiii ,k+q v /2.n2,k-q v /'2 
n '\ >k' — q y /2,n' 2 ,k' +q y /2 



(7.1) 



for the conventional quadratic electron spectrum (since 
only the e = LL states are involved). Here, 
rii, ri2, rii, n' 2 are the LL indices, V(q) is the interaction 
potential in the momentum space, q = |q|, and 



K nn t(q) = 



+ 00 



dxe" 



are the form- factors, with 4> n {x) the magnetic oscillator 
eigenstates. 



In Ec, <2.3), Vg£ = and 



V{q) 



V Q (q) 



i + n(?)v (g) 



is the screened Coulomb potential, with Vo(q) 
27re 2 / (nq) and H(q) the static polarization operator. 



In Eqs. (|2.4D, (|6.3|>, and <6.4|, Vgg = V$$% 



and 

V(q) = 1 is the point interaction potential of unit 
strength. 
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